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DEFLECTION  AND  STRESS  ANALYSIS  OF  THIN  SOLID  WINGS  OF  ARBITRARY  PLAN  FORM 
WITH  PARTICULAR  REFERENCE  TO  DELTA  WINGS  1 

1  By  Manuel  Stein,  J,  Edward  Anderson,  ami  John  M.  Hedgepeth 


SUMMARY 

The  structural  analysis  oj  arbitrary  solid  cantilever  wings  by 
small-deflection  thin-plate  theory  is  reduced  to  the  solution  oj 
linear  ordinary  differential  equations  by  the  assumption  that 
'  the  chordwise  deflections  at  any  spanwise  station  may  be  ex¬ 

pressed  in  the  form  oj  a  power  series  in  which  the  coefficients  are 
Junctions  oj  the  sjranwise  coordinate.  If  the  series  is  limited  to 
the  first  two  and  three  terms  (that  is,  if  linear  and  parabolic 
chordwise  deflections,  respectively,  are  assumed >,  the  differential 
equations  for  the  coefficients  are  solved  exactly  for  uniformly 
loaded  solid  delta  wings  of  constant  thickness  and  oj  symmetrical 
double-wedge  airfoil  section  with  constant  thickness  ratio.  For 
cases  for  which  exact  solutions  to  the  differential  equations 
cannot  be  obtained,  a  numerical  procedure  is  derived.  Experi¬ 
mental  deflection  and  stress  data  for  constant-thickness  delta- 
plate  specimens  of  1,5*  and  G0°  sweep  are  presented  and  art 
found  to  compare  favorably  with  the  present  theory. 

INTRODUCTION 

One  of  tho  present  trends  in  the  development  of  high-speed 
airplanes  and  missiles  is  toward  tho  use  of  thin  low-aspcct- 
ratio  wings.  Tho  structural  analysis  of  these  wings  often 
cannot  he  based  on  beam  theory  sinco  tho  structural  defor¬ 
mations  may  vary  considerably  from  those  of  a  beam  and, 
indeed,  may  more  closely  approach  thoso  of  a  plate.  In 
cases  where  the  wing  construction  is  solid  or  nearly  solid  tho 
uso  of  [date  theory  in  tho  analysis  is  particularly  valid,  and 
it  is  this  type  of  wing  which  is  considered  in  the  present 
report. 

j  Exact  solutions  to  tho  partial-dilfcrcntinl  equation  of  plate 

theory  arc  not  readily  obtained,  especially  for  plates  of 
1  arbitrary  shape  and  loading;  however,  a  number  of  approxi- 
i  mate  solutions  to  specific  problems  on  cantilever  plates  havo 

appeared  in  tho  literature  (see,  for  example,  refs.  1  to  7). 
Of  the  approaches  used  in  these  references,  only  the  one  in 
references  6  and  7  is  readily,  applicable  to  plates  of  arbitrary 
plan  form,  thickness  distribution,  and  load  distribution;  thus 
it  is  the  most  useful  one  for  the  analysis  of  actual  wings. 

In  reference  0  the  cantilever-plate  problem  is  simplified  by 
the  assumption- that  tiic  deformations  of  tiio  plate  in  iiio 
chordwise  direction  (parallel  to  the  root)  are  linear.  By 
minimizing  tho  potential  energy  of  tho  plate,  tho  partial- 
differential  equation  of  plate  theory  is  replaced  by  two 


simultaneous  ordinary  differential  equations  for  the  spanwise 
variations  of  the  bonding  deflection  and  twist.  In  reference  7 
the  same  ordinary  differential  equations  arc  obtained  in  a 
different  manner.  Refinement  of  the  analysis  by  inclusion 
of  the  effect  of  parabolic,  cubic,  or  higher-order  chordwise 
camber  terms  is  indicated  in  reference  0,  and  ns  the  order,  of 
refinement  is  increased  a  corresponding  increase  in  the  num¬ 
ber  of  ordinary  differential  equations  is  obtained. 

In  tho  present  report,  which  is  an  extension  of  reference  0, 
a  general  set  of  ordinary  differential  equations  is  presented 
which  may  bo  used  to  obtain  any  desired  degree  of  approxi¬ 
mation  to  tho  deflection  of  the  plr.tc.  These  equations  are 
solved  exactly  for  several  cases  of  delta  plates  under  uniform 
load  first  by  considering  linear  chordwise  deformation  only 
and  second  by  including  tho  effect  of  parabolic  chordwise 
camber.  Comparisons  are  drawn  between  the  stresses  and 
deflections  computed  from  tho  equations  of  each  approxi¬ 
mation  and  ulso  with  some  experimental  results. 

The  differential  equations  presented  contain  coefficients 
that  depend  on  the  plan  form  and  stiffness  distribution  of  the 
plate  and  on  the  loading.  In  this  report,  the  plates  con¬ 
sidered  in  detail  linvc  coefficients  such  that  tho  differential 
equations  enn  be  solved  exactly ,  however,  in  cases  for  which 
exact  solutions  cannot  be  obtained  a  numerical  procedure 
must  bo  used.  One  such  procedure  is  derived  and  its 
accuracy  is  demonstrated. 

SYMBOLS 

l  length  of  pinto  measured  perpendicular  to  root 

c  root,  chord  of  pinto 

p  lateral  load  per  unit  area,  positive  in  s- 

dircction 

l  local  thickness  of  plate 

i„  nverago  thickness  of  plate 

El1 

l)  local  flexural  stifTncss, 

U  flexural  stiffness  based  on  average  thickness, 

EU 

12(1 -u?) 

E  modulus  of  elasticity  of  material 

fi  Poisson’s  ratio 

w  deflection  of  plate,  positive  in  z-dircction 

x,  y,  z  coordinates  defined  in  figure  1 


*  Supersedes  NACA  TN  262),  "Deflection  and  Stress  Analysis  of  TliLn  Solid  Wings  ol  Arbitrary  Plan  Form  With  I'articukir  Reference  to  Delta  Wings"  by  Manuel  Stein,  J.  KdwrurU 
Anderson,  and  John  M.  Hedgepeth,  IM2, 
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Kiouke  1.— Coordinate  system  used  in  tlio  present,  analysis  fora  canti¬ 
lever  plate  of  arbitrary  sliape  with  arbitrary  thickness  variation. 


to 

function  of  x,  coefficient  in  power  series  for 

.V 

deflection  w—'P.  <o.f z)>/* 

»-o 

ci  Mi  cjM 

functions  defining  plan  form  (sec  fig.  1 J 

Xl 

variable  obtained  by  transformation  7,=1 — j 

CXf  <Jy 

normal  stresses 

TZV 

shear  stress 

a 

maximum  principal  stress 

X 

l  /3 

aspect-ratio  parameter,  (1  —  p) 

RESULTS 

The  derivation  of  the  general  set-  of  ordinary  differential 
equations  is  given  in  appendix  A.  The  general  procedure 
outlined  in  reference  G  is  followed;  that,  is,  the  deduction 
of  the  plate  w  is  expanded  into  a  power  series  in  y  the  chord- 
wise  coordinate  with  coedicicnts  which  are  functions  of  x 
the  spanwise  coordinate  (see  fig.  1) 

w = v>oM + pi(*)i/+ *>»(*)»*+  .  .  .  -rv>.v(x)!/-v  (1) 

Equation  (1)  is  substituted  into  the  expression  for  the  po¬ 
tential  energy  of  the  plate-load  combination  which  is  in 
turn  minimized  by  the  calculus  of  variations  with  respect 
to  each  of  the  coefficients  *>„.  The  results  of  the  variational 
procedure  appear  as  Ar-}-l  simultaneous  differential  equa¬ 
tions  with  the  coefficients  v>«  as  unknowns. 

By  taking  a  sufficient  number  of  terms  in  tho  expansion 
of  w,  the  resulting  differential  equations  can  be  used  to  ob¬ 
tain  any  desired  degree  of  accuracy  in  the  solution  for  the 
dcdections  of  any  given  cantilever  plate  subjected  to  an 
arbitrary  lateral  load.  Of  most  interest,  perhaps,  are  the 
particular  cases  for  N=  1  and  N=2,  which  arc  obtained 
from  the  general  set  of  equations  and  arc  simplified  in 
appendix  A.  The  case  for  Ar=i  (also  derived  in  refs.  0 


and  7)  includes  linear  ehordwise  deflections,  and  the  caso  for 
A'=2  takes  into  account  parabolic  ehordwise  curvature. 
Although  for  most  practical  problems  tho  solution  by  the 
parabolic  theory  should  be  adequate,  eases  might  exist  in 
which  cubic,  quartie,  or  even  higher-order  chord wiso  terms 
should  bo  included,  depending  on  tho  convergence  of  tho 
series  for  tho  configuration  considered. 

The  particular  equations  for  A7— l -and  Ar=2  are  used  to 
determine  the  deflections  and  stresses  of  tho  following  can¬ 
tilever  plates  subjected  to  uniform  lateral  load: 

(1)  A  45°  delta  plate  of  uniform  thickness 

(2)  A  GO0  delta  plate  of  uniform  thickness 

(3)  A  45°  delta  plate  of  symmetrical  double-wedge  airfoil 
section  with  constant,  thickness  ratio 

Fortunately,  for  these  configurations,  the  solution  can  be 
carried  out  exactly  by  both  the  linear  and  parabolic  theories, 
and  the  details  cf  these  exact  solutions  are  included  in 
appendix  B.  In  general,  however,  exact  solutions  cannot 
bo  obtained  and  some  numerical  method  must  be  used. 
One  such  method,  based  on  replacing  derivatives  by  their 
first-order-approximation  difference  forms,  is  derived  in 
appendixC. 

A  summary  of  the  results  for  the  three  particular  problems 
is  shown  in  figures  2  to  11.  Deflections  obtained  by  the 


0  .2  .4  .6  8  10 
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Fiourk  2. — Deflections  of  a  45°  delta  plate  of  uniform  thickness  1 
uniform  load.  ' 
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linear  theory  and  the  parabolic  theory  for  the  three  con¬ 
figurations  arc  compared  in  figures  2,  3,  and  4.  Stresses 
obtained  by  the  linear  theory  and  the  parabolic  theory 
for  the  three  configurations  aro  compared  in  figures  5,  0, 
and  7.  Where  available,  experimental  deflections  and 
stresses  arc  also  shown  in  these  figures.  The  details  of  tho 
procedure  used  to  obtain  tho  experimental  deflections  of 
tho  45°  and  00°  uniform-thickness  plates  and  tho  experi¬ 
mental  stresses  in  the  45°  uniform-thickness  plate  arc  con¬ 
tained  in  appendix  D,  whereas  the  experimental  root  stresses 
for  the  00°  uniform-thickness  plate  were  obtained  from 
reference  8.  figures  8  to  11  present  the  comparison  be¬ 
tween  deflections  and  stresses  computed  from  the  exact 
solutions  of  tho  differential  equations- and  those  computed 
from  the  numerical  solutions  of  the  same  equations. 

DISCUSSION 

Tho  results  shown  in  figures  2  and  3  indicate  that,  with 
regard  to  deflections,  cither  the  linear  theory  or  the  parabolic 
theory  is  adequate  for  the  case  of  a  constant-thickness  delta 
plato  subjected  to  a  uniform  load,  the  comparison  being  some¬ 
what  better  for  tho  00°  plate  than  for  tho  4.r>°  plate.  If 
accurate  slopes-in  the  chordwiso  direction  (angle  of  attack) 


o  z  .4  .6  .0  10 

y/c 

Figure  3. — Deflections  of  a  G0°  delta  plate  of  uniform  thickness  under 
uniform  load. 


are  desired,  how  ever,  the  parabolic  theory  must  be  used 
because  the  error  in  the  angle  of  attack  as  computed  by  the 
linear  theory  is  ns  much  ns  30  percent  (see  figs.  2  and  3). 
The  appreciable  nnticlastic  curvature,  evidenced  by  the 
cxpct  .mental  results  of  figures  2  and  3,  may  be  important 
aerody  nanueally  and  is,  of  course,  nut  taken  into  account  by 
the  linear  theory. 

The  apparent  convergence  of  the  aforementioned  series 
in  the  case  of  the  double-vv  edge-section  pinto  (see  fig.  4) 
implies  that  the  linear  theory  is  adequate  foi  this  case.  The 
lack  of  chordwiso  curvature  in  tho  result  obtained  by  the 
parabolic  theory  is  attributable  to  the  fact  that  tho  natural 
tendency  of  the  plate  to  have  nnticlastic  curvature  is  canceled 
by  the  opposite  tendency  of  tho  thin  edges  to  bend  down 
tinder  the  load.  Unfortunately,  no  experimental  results 
are  available  for  (his  configuration. 

Tn  figure  4  tho  plate  stiffness  U  in  tho  nondimensionnl 
parameter  wlj/pl'  is  the  local  value  of  1)  at  a  point' where 
the  thickness  is  equal  to  the  average  thickness  of- the  plato 
as  a  whole.  Thus  the  results  of  figure  4  are  comparable 
with  the  results  of  figure  2  on  an  equal-weight  basis.  V 
can  lie  seen  that  the  deflections  of  the  doublc-wcdge-section, 
constant-thickness-ratio  plate  arc  everywhere  less  than 


Figure  4. — Deflections  of  a  45°  delta  plato  of  symmetrical  dotilde- 
ivcdgoairfoil  section  and  constant  thickness  ratio  underuniform  load. 
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those  of  tho  uniform-thickness  pinto  Although  they  increase 
rapidly  near  the  tip.  This  curling-up  or  singularity  in  slope 
at  the  til>  is  a  result  of  using  a  smnll-dcllection  theory  and 
probah'y  would  not  he  so  marked  in  an  actual  case. 

Tho  stress  results  for  tho  45°  and  00°  uniform-thickness 
delta  plates  indicate  that  both  the  linear  and  tho  parabolic 
theories  arc  adequate  and  that  tho  parabolic  theory  is  bettor 
than  tho  linear  theory  only  near  tho  root.  It  should  lie  noted 
that,  although  the  maximum  principal  stress  over  a  large 
part  of  tho  45°  plate  is  plotted  in  figure  5,  only  the  stresses 

normal  to  tho  root  along  the  lino  j=0.0087  of  tho  00°  plate 

are  plotted  in  figure  0  since  only  those  stresses  are  giren  m 
reference  S.  (The  maximum  principal  stress  and  the  stress 
normal  to  the  root  are  theoretically  equal  at  the  root  since 
the  root-  shear  stress  is  zero.) 

Experimental  data  are  lacking  for  the  doublc-wedge-section 
delta  plate  and,  therefore,  only  theoretical  stresses  are  shown 
in  figure  7.  As  in  the  case  of  deflections,  the  results  obtained 
from  the  linear  theory  and  those  obtained  from  the  parabolic 
theory  are  almost  coincident,  the- difference  being  greatest 
near  the  root.  Figure  7  has  also  been  plotted  so  that  tho 
results  are  directly  comparable  with  those  for  the  45° 


I'Tuuin:  5. — Maximum  principal  stresses  in  n  15°  delta  plate  of  uniform 
thickness  under  uniform  load. 


uniform-thickness  plate  in  figure  5  on  an  c  unl-weight 
basis.  As  can  be  expected,  the  doublc-wedge-section, 
constant-thichness-ratio  plate  is  a  better  design  stiucturally  , 
the  sties.  »s  in  the  donble-wedge-scction  plate  are  everywhere 
smaller  and  are  almost  constant  in  the  spnnwiso  direction. 

The  theoretical  results  in  figures  2  to  7  have  been  obtained 
from  exact  solutions  of  the  differential  equations  of  the 
linear  and  parabolic  theories.  In  order  to  test  the  reliability 
of  the  numerical  method  derived  in  appendix  C,  the  differ¬ 
ential  equations  were  also  solved  numerically.  The  results 
shown  in  figures  8  and  9  indicate  that  the  agreement  is  good 
between  tho  numerical  solution  in  which  five  equal  intervals 
were  used  ami  the  exact  solution  of  the  differential  equations 
for  the  enso  of  the  45°  uniform-thickness  plate.  The  same 
good  agreement  can  he  expected  in  other  cases  where  the 
thickness  and  load  distributions  arc  not  too  erratic  and  where 
the  plntc  stiffness  does  not  go  to  zero  at  the  tip  that  is, 
when  no  singularities  appear  at  the  tip. 


n  GO’  delta  plate  of  uniform  thickness  under  uniform  load. 
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Sinco  the  efficacy  of  tlio  numerical  method  depends  on 
liow  well  parabolic  arcs  fit  tlio  various  functions  between 
stations,  serious  error  can  result  from  blind  application.  An 
examplo  of  the  seriousness  of  these  errors  and  of  tho  manner 
in  which  they  can  bo  remedied  is  shown  in  figures  10  and  11. 
In  these  figures  a  comparison  is  made  between  exact  and 
numerical  results  obtained  on  tho  45°  doublc-wedgc-scction, 
cpnstant-thickness-ratio  plate.  As  can  be  expected,  the 
fire-station  numerical  solution  fails  to  follow  the  oxact  solu¬ 
tion  in  tho  neighborhood  of  the  singularity  at  tho  tip.  Sinco 
tho  region  of  trouble  is  localized  at  the  tip,  a  reasonable 
remedy  would  be  to  dccreaso  tho  spacing  of  the  station  points 
near  tho  tip.  This  dccreaso  in  spacing  may  be  accomplished 
cither  by  using  a  greater  number  of  equally  spaced  stations 
or  by  using  unequally  spaced  stations  crowded  near  the  tip. 
The  increase  in  accuracy  obtained  by  increasing  the  number 
of  equally  spaced  station  points  to  ten  is  shown  in  figures 
10  and  11. 

CONCLUDING  REMARKS 

The  general  method  presented  herein  for  finding  deflec¬ 
tions  and  stresses  of  solid  or  nearly  solid  wings  is,  in  principle, 


Fluent.  7. — Maximum  principal  stress  in  a  45°  della  plate  uf  e.i  iniuct- 
ncal  doiiltlo-wcdge  airfoil  section  and  constant  thickness  ratio  under 
uniform  load. 
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capable  of  yielding  arbitrarily  accurate  results  for  any  con¬ 
figuration.  It  is  seen  thnt,  for  the  examples  considered, 
only  the  first  two  or  three  terms  m  the  series  expansion  need 
be  considered  to  obtain  adequate  accuracy.  In  addition, 
for  most  practical  plate-like  wings  with  clnmpcd  roots  the 
first  two  or  three  terms  will  probably  bo  adequate,  although 
problems  may  exist  wherein  moro  terms  are  needed. 

The  numerical  procedure,  derived  for  application  in  cases 
whero  exact  solutions  cannot  be  obtained,  gives  good  agree¬ 
ment  when  compared  w  ith  exact  solutions  if  enough  stations 
aro  taken  along  the  spnn.  The  necessary  number  of  stations 
is  dependent  on  the  type  of  thickness  nml  loading  distribution 
considered,  five  equally  spaced  stations  being  enough  for 
the  uniform-thickness  delta  wing  subjected  to  uniform 
loading  and  ten  being  necessary  for  tho  double-wedgc-scction, 
constant-thickness-ratio  delta  wing  subjected  to  uniform 
loading. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va,,  November  SO,  1051. 


Fluent  8.  Numerical  ami  c.xacl  solutions  uf  the  differential  equations 
fur  the  deflections  uf  the  free  edges  of  a  45°  delta  plate  of  uniform 
thickness  under  uniform  load. 
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I'loURB  11. — Numerical  and  exact  solutions  of  the  differential  equa¬ 
tions  (obtained  by  assuming  linear  chordniso  deflections)  for  the 

maximum  principal  stress  along  the  line  *>f  ft  In0  delta 

plato  of  symmetrical  double-wedge  airfoil  section  ami  constant 
thickness  ratio  under  uniform  load. 


APPENDIX  A 


DERIVATION  OF  DIFFERENTIAL  EQUATIONS 


Tito  structure  considered  herein  is  n  thin,  elastic,  isotropic, 
cantilever  plate  of  arbitrary  plan  form  and  slowly  varying 
thickness  subjected  to  distributed  lateral  load  (see  fig.  1). 
By  assuming  that  the  deflection  of  the  plato  can- be  repre¬ 
sented  by  a  power  series  in  the  chordwisc  coordinate  and 
by  applying  tho  miniinum-potcntial-energy  principle,  aset 
of  ordinary  differential  equations  in  the  spamvisc  coordinate 
is  obtained  from  which  tho  coefficients  of  the  power  series 
may  be  determined.  First  the  general  sot  of  equations  is 
derived;  then  tho  particular  equations  for  tho  eases  of  linear 
chordwisc  deflections  and  parabolic  chordwiso  deflections  are 
deduced  from  the  general  set  and  simplified. 


General  equations.  Tho  potential  energy  of  the  system 
under  consideration  is 


. . 

2*  5  §*+ 2(1  - ^ </y dx  (A1) 


in  which 


D{x,y) 


12(1  -ri 


and  p(x,y)  is  tho  distributed  lateral  load. 


M3M-U — } 
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The  assumption  is  made  that,  (ho  dolicctiou  to  can  he  represented  by  the  power  series 

w=S  *>.(*)</■  (A2) 

tl-U 

Substitution  of  this  expression  for  to  into  cquntion  (Al)  gives 

J'l  ( J  x  x 

0  \ “  ro-0  n«u 

2pn(n— l)«»+.-,p,»'V,-r2(l— M)wnaB+1,_,p1.V.'|— 2jI».+iV>.|  (A3) 

n-u  ) 

in  which 

«,-f  fJ'>  (r-  1,2,  .  .  .  2A'+1) 

7>’  JeM  p{*,y)>f~'ily  (r=  1,2,  .  .  .  A'+l) 

and  tlio  primes  denote  differentiation  with  respect  to  x. 

Minimi/', ation  of  the  potential  energy  l>y  means  of  the  calculus  of  variations  gives 

{(Potential  energy) =0 

=  f  dx  S  £  [«»+.4i(v»i"*V«"+¥>«"iv»™")+w,«("i  —  1)  (»—  l)<i*+»-)(s»m5v’»+v’,{p«)Hh 
Jo  \  “  n»“0  n»0 

2jiti(tt— l)«M+,.1((s„,'{^114-v»J^r,")+2(l— n)mn«l,4..,(v>«'^»,+«>i1'{p,,')l-’2j7'»+i^«| 

n-e  ) 

Integrating  by  parts  and  collecting  terms  results  in 

0=  (  </*£{*>„  {s  I(«*+.+iv>™")"+ft»'(tn-l)  (<i»+,-iv>»)"-2(l-it)m/t(«w+,.,v>n')'+><«(«-l)«»+»-iv»»i"+ 

Jo  n-0 

mn(m— 1) (a -1) «„+«., e>J-7t,+,|  +  {Siy.'  S  (<'«+.+iv>n"+fi»i("i-l)«i™+.-ii>»ll  -{zifa.  £  I(«m+»+iv>«")'+ 

)  VH“U  W-0  JO  v«“0  M«0 

fit/i  (m  — 1)  y„)'—2  (l-^)ot  h<i„+,. i  •  (AO) 

Everywhere  in  the  region  of  the  plate,  except  at  the  boundary  r=0,  the  variation  of  to  is  arbitrary.  At  *=()  the 
cantilever  boundary  conditions 

Dto  . 

*-ss-° 

yield 

<A.(0)-p„'(0)=0  (n=0, 1, ...  A')  (AO) 

and  therefore  the  variation  in  these  quantities  must  also  bo  zero  at  /=(). 

Equation  (AO)  is  then  satisfied  if,  in  nddition  to  equation  (AO), 

2l(fl«+»+iP«i*)"+t‘m(",—  l)(«-»*»-i¥>r»)"— ^ 2(1—  ii)mn(a„i.n.lvmy+.iin(n— 1)«„+._,*>„"+ 

fn*0 

m/i(m—  l)(;t-I)a„+,_3v>„]=7t,+l  (n=0, 1, . . .  A')  (A7) 

Sl"».+.+iv’«"+n'«('«— l)ar,+,-wJz.i-0  (»= 0. 1,  ...AO  (A8) 

W"0 

and 

S((««.s»+a»«")'+n«('n-l)  2(1  -^)w««„+,_, =0  (n= 0, 1, . . .  A’)  (AD) 

w»o 

Equations  (A7)  form  a  set  of  Ar|  1  simultaneous  ordinary  differential  equations  for  the  functions  <fn(x).  The  functions  <f, 
are  completely  determined  by  these  differential  equations  and  the  boundary  conditions  (AO),  (A8),  and  (A9). 
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DEFLECTION  AND  STRESS  ANALYSIS  OF  THIN  SOLID  WINGS  OF  ARBITRARY  PLAN  FORM  0 


Particular  cate  of  linear  chordwise  deflections. — If 
A'=  1 ,  the  deflection  function  becomes 

u>=Vo+?/fi  (A10) 

n  linear  function  in  the  chordwise  direction,  where  vo  is  the 
bending  deflection  and  vi  is  the  twist.  Equations  (A7) 
become 

(«.«*")"+  («.*.")"=/»  (All) 

(A  12) 

Tiie  root  boundary  conditions,  given  by  equation  (AO), 
becomo 

ft(0)=vo'(0)==v.(0)==v>.'(0)=0  (A  lit) 

'I'lio  tip  boundary  conditions,  given  by  equations  (AS)  and 
(A9),  becomo 

(<inV'+<hv'i")x-i=0  (AH) 

(a4v\,"-(-(ij(5,"),-i=0  (Aid) 

[(OivV')'+(«1*»i"),).-.=0  (A10) 

l(aJvV')'+(oav*l")'-2(l-M)<ilv.1'),.i=0  (A  17) 

Equations  (All)  to  (A  17)  arc  the  differentia!  equations  and 
corresponding  boundary  conditions  presented  in  reference 
0  (if  only  distributed  load  is  considered)  where  the  symbols 
11' and  6  arc  used  instead  of  vo  and  vi,  respectively. 

If  equation  (All)  is  integrated  twice  and  the  boundary 
conditions  (AH)  and  (AlG)  nro  used, 

<*"=-5>"4“  ■j'j'r.M  (A  is) 

Substitution  of  po"  into  equations  (A  12),  ami  (A  17) 

gives 

(6,^")"-2(1-p)(«1^')'=;>1-(2|J/[V, </*’)"  (Alt)) 

(4h9i"),.i=0  (A20) 

l(W')'-2(l -,/)«, *>,%.,==<>  (A21) 

in  which 


If  equation  (A  19)  is  integrated  once  and  tho  boundary 
condition  (A21)  is  used, 

(4iv5i")'-2(l~/i)(iiy>i'=  — J  />,  dz2)  (A 2 2) 

The  differential-equation  (A22)  is  a  second-order  differential 
equation  in  vi'-  The  twist  vi  and  then  the  bending  deflection 
Vi)  aro  obtained  by  solving  equations  (A22)  and  (A18),  respec¬ 
tively,  by  applying  tho  boundary  conditions  (A13)  and  (A20). 

Particular  case  of  parabolic  chordwise  deflections.— The 
effect  of  parabolic  chordwise  camber  may  bo  included  by 


letting rV=2  in  the  general  powerseries  (eq.  (A2)).  If  ,V=2, 
the  deflection  function  becomes 

I/!--  Vo4->/vrK'/V2 

Hero  ift  represents  the  spnnwiso  distribution  of  parabolic 
chordwise  camber.  For  this  case  the  differential  equations 
(A7)  become 

(«iVo”)"  4-  (<ijVi”)*+("ji9i,’)*+2a('iiVi)"  ~Pi  (A23) 

2(l-#.)l(«iv-,')'+2(«lV,J')'l=^  (A24) 

(««*")"  -i-(«iv>i*)" + («4pi")'' +2r[«i  w'  d-rtivi*  +«jv»i"  4- 

(ojVi)*)- 4(1 — ri)[(<'2ifI(,),+2((/jv/),l+-lfliVj  ~ /'i  (A25) 

with  the  boundary  conditions 

*>(0)  =w>'(0)  =  if  I  (0)  —V,'  (0)  =  Vi(0)  =¥>/«>)  -=0  (A2t>) 

(,,lVu"‘r'liVl”+('lV:”4‘2p'liV‘))/-l— b  (A27) 

(«rv>o"+<i>t>i'',4*«iv''j"-{'2jiaivi),.|— 0  (A28) 

('iivo"+«iv'i"+"sVi"-t-2V«jVi))-i=0  (A29) 

l(<ii^i),,),4-(«2v!i”),+ (o»^/)'+2n(«iV})'),.i=0  (A:iO) 

l(«2Ya") '  +  («av»i") ' + 4- ' — 

2(1 -*)(«, *,'+2cW)|,.l=0  (AH  1) 

[(«3V-o")'+ («.*,")'4- (wSY+MomY- 

4(1— h)(h2vi'+2hjvV)1i-i=:0  (A;i2) 

If  equation  (A23)  is  integrated  twice  and  the  boundary 
conditions  (A27)  and  (A30)  are  used, 

f  f  Pi  (lx*  (A33) 

”1  ”1  J I  J  l 

Substitution  of  vs"  into  the  remaining  differential  equations 
and  boundary  conditions  results  in 

(&,V.T+(W)"-2(1  -p)  |(«.v,')'4-2(<W)'l 

(A34) 

(btvi)"  4  ifri'pi)"— 4(1—  /»)  [(«sv>iO' 4  2(«jV2,),]4- 4(1  —  /i*)«iV2 

=,'3“2" XT,,<  w-(U/,Tft  -y  (A35) 


(('iv>i"4-6jv’j")i.i=0  (A3C) 

1(&i¥,iV+  )'  —  2(1  —  p)t<l|V>i*  +  2«2V2'>li.i— 0  (A37) 

(biipi"+bypi")i^i=0  (A3S) 

Wifi")'  — 4(1  —  i2)((i2^i,4-2«3V2,))i-i=0  (A39) 

Vi  (0)  — -  V! '  (0)  —  v>i(0)  =V2'(0)  — 0  (.440) 
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Thus  ci  and  ct  arc  obtained  by  solving  equations  (A35)  and 
(A41)  with  tho  boundary  conditions  (A36),  (A38),  (A39),  and 
(A40).  Subsequently,  c»  can  be  obtained  by  solving  equa¬ 
tion  (A33)  witli  tho  boundary  conditions  co(0)=v>j'(0)=0. 


Stresses. — After  the  approximate  deflection  of  tho  plate 
is  determined  from  equations  (Al8)  nnd  (A22)  or  from  equa¬ 
tions  (A33),  (A35),  and  (A41),  the  extreme-fiber  stresses  may 
bo  calculated  from  the  well-known  equations  of  thin-platc 
theory,  which  are  (see,  for  example,  ref.  9): 

C  D  /b’w  .  bho\ 

a'~ir\Tx*+,1w) 

6  D  /bhu  ,  bho\ 

_G(1— ix)D  bho 
r"~  t‘  bxby 

The  maximum  principal  stress  <r  at  any  point  in  the  plate 
can  bo  determined  from 

O  =  ±  £  V(«Tj — <i»)M-4  r„2 


APPENDIX  B 


EXACT  SOLUTIONS  OF  DIFFERENTIAL  EQUATIONS  FOR  SOME  SPECIFIC  DELTA-PLATE  PROBLEMS 


Tlio  differential  equations  of  appendix  A  for  linear  and 
parabolic  chordwise  deflections  arc  solved  exactly  for  uni¬ 
formly  loaded  delta  plates  of  constant,  thickness  and  of  sym¬ 
metrical  double-wedge  airfoil  section  with  constant  thickness 
ratio.  The  equations  for  deflections  obtained  by  tho  linear 
theory  are  presented  in  terms  of  tho  aspect-ratio  parameter 
X  for  both  kindtf  of  delta  plates.  The  equations  for  deflec¬ 
tions  obtained  by  the  parabolic  theory  are  presented  for 
l  1 

-=  1  and  — -  with  for  the  constant-thickness  delta  plate 

Co  o 

and  for  -=  1,  also  with  for  the  delta  plate  of  symmet- 
C  o 

rical  double-wcdgo  airfoil  section  -with  constant  thickness 
ratio. 

If  the  x-nxis  is  passed  through  the  edge  perpendicular  to 

j 

the  root  and  tho  substitution  *i=l— y  is  made,  the  differen¬ 
tial  equations  arc  clearly  of  the  homogeneous  type  for  which 
the  solutions  arc  of  the  form  xA,  where  y  is  a  constant.  For 
the  configurations  considered,  tho  functions  that,  define  the 
plan  form  (sco  fig.  1)  are  then  ci(z)=0  and  ci(x)=cxh  where  c 
is  the  root  chord.  In  all  tho  equations  of  this  appendix  tho 
primes  denote  differentiation  with  respect  to  tho  now  inde¬ 
pendent  variable  rt. 

DELTA  PLATE  OP  UNIFORM  THICKNESS  UNDER  UNIFORM  LOAD 

Since  the  stiffness  D  is  a  constant  for  uniform-thickness 
plates,  the  coefficients  in  tho  differential  equations  (seo  cq. 
(A4))  become 

a«=~T“  (Bln) 


chordwise  deflections  result: 


tovvy-ioxw*— 2-j 

H  e  „  .  pl{  , 
<Po  —  +6  dx* 


k=,c  Vf  (1~") 


Tho  boundary  conditions  to  he  used  with  tlicso  equations  are 
obtained  from  equations  (A13)  and  (A20)  and  are 

v»(l)=vVtl>=v.(l)=v>.'(l)=0  (B4) 

(J'.Vi"),1-o=0  (B5) 

Tho  general  solution  of  equation  (B2)  is 


Vi  =Aixi'~t+AtX\~'~'- 


*1*  pi' 
"4(1 -2X0  Dc 


Y=Vr+ 1 GX2 

and  -l,  and  A2  are  arbitrary  constants.  Since  X*  is  inherently 
positive,  the  boundary  condition  (1)5)  requires  that  As=0. 
One  integration  of  equation  (1)0)  and  the  application  of  the 
conditions  y>i(l)~ s^i'O)— 0  yields 


4(1—  2X*j  Dc 


a«— x<"  (Bln)  if  equation  (1)3)  is  solved  for  y"o  with  tho  conditions 

v,»(l)=vt,(l)=0.  the  result  is 

5,=a,-2A=0|  *,»  (Bib)  pl»  1  T2  f  1 

a,  12  *»-g2Ji‘=33?Lol°”4X,kl "T~)~ 

6J==u,_2i2l=®|Il.  (Bic)  (B8) 

(Bpi)  Substitution  of  equations  (1)7)  and  (1)8)  into  the  equation 

w=ifl,+t/ip, 

_ pcn  „  -  . 

V*~  ~  *i  (Die)  gives  the  expression  for  the  deflection  v>  of  tho  plate  under 

tho  assumption  of  linear  chordwise  deflections. 

Solution  for  linear  chordwise  deflections,— If  the  co-  Solution  for  parabolic  chordwise  deflections— If  the  co¬ 
efficients  given  by  equations  (Bl)  are  substituted  into  cqua-  efficients  given  by  equations  (Bl)  are  substituted  into  equa¬ 
tions  (A22)  and  (AIS)  and  the  independent  variable  is  (ions  (Ally,  (A35y,  and  (A.53)  and  the  independent  variablu 

changed  to  the  following  equations  for  linear  is  again  changed  to  j-|  =  I  —  j>  the  following  equations  for 


Solution  for  linear  chordwise  deflections, — If  tho  co¬ 
efficients  given  by  equations  (Bl)  are  substituted  into  equa¬ 
tions  (A22)  and  (A18)  and  the  independent  variable  is 


i 
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parabolic  chordwise  deflections  result:  |  The  general  solution  is  the  sui 


The  general  solution  is  the  sum  of  the  homogeneous  solu- 

(nW+faWr-  icx‘(W+*,W>=-2  g  «,*  ,!o"s  amI  ,h0  ',arlic,,l“r  inU‘gral 

foWT'+Jf  (z.wr- 1  ex*  [(*,V)'+§  (X,’e^')']+  7 

,  ,  (1317) 

04..  H-ji  4/_,  2  n\}\pl{  , 

-g-X‘ x,Cfj=-Tj (<  (Bio) 

,  ,  where  the  values  7,  are  the  six  roots  of  the  characteristic 

vfo"=—  £ i\<pi" — -■  xiVj"_ Sp/Vi+Tj- rr  (Bll)  equation  (131G)  and  the  coeflicients  /l.  and  /3„  are  the  co- 

'  J  1  efficients  corresponding  to  each  of  these  roots.  After  intc- 

Tho  boundary  conditions  to  bo  used  with  these  equations  arc  gration  <ci  becomes 

*(l)-w/O)“«(l)-*'(l)-«<l)-w'(l>*=0  (1312)  ¥>i=2/1„  ^+xl,4^+/l(  03JS) 

n*»l  In  O 

(^W/4-^lV1',)x.-0=0  (13 1 3)  11*’/  /  *«  /n.,x*  /  1 

I  ho  general  solution  for  ^  from  equation  (Bll)  is  found 

(*.Vi"+77  =0  (1314)  tobc  6 

v  y>«*  2  C^- * ' + CPxt'+C,x,  -I-  0,  (1319) 

[~(yiVi")'+T7  toWT-iox'fx.V+l  =0 

L  \  o  /Jn-o  wlieic,  for  »=  1, 2, .  .  .  0, 

(1315) 

The  homogeneous  solutions  of  the  simultaneous  equations  C„=—  r;.-,  {-tt  *  *l»+f|(Ti.—  1)(t«— -)+|,^r"l^«| 

(139)  and  (1310)  arc  of  the  form  M  '  L  ' 


n«l 


n**l 


v.=2*l,ir-+^e;r+/l#  03JS) 

n-l  T11  o 

The  general  solution  for  y\>  from  equation  (1311)  is  found 
to  be 

^=2C,x,».'"-1-(7,x1*+C5x1-1-C',  (1319) 

11-1 

whet c, -for  n=l,  2, .  .  .  fi, 


V>|'=AX|T“‘ 

ipi—lixt'-' 


O- h['1>+?(,+i-;>-««J 


Substitution  of  these  expressions  into  the  homogeneous  parts  The  coeflicients  At  to  A,„  A„  C„  and  C,  must  he  determined 

of  equations  (139)  and  (1310)  leads  to  the  following  chnrac-  by  the  boundary  conditions  (1312)  to  (1315). 

tcristic  equation  from  which  X  may  be  determined:  -A  complete  set.  of  coefficients  is  given  in  the  following  table 

[for  delta  plates  with  Poisson’s  ratio  n  equal  to  1/3  and  with 

320  (4+  'tA‘)X«+480XM-0 It’—  ,  l  ,  .  ^3  „  „  ,  „  .  , 

\  1—^/  J  X=-=l  and  Deflection  curves  plotted  from  these* 

4  [l280  {±£  X‘+S0  (4+1+7)  X*+9GX2+ 1]=0  ,CSH"S  ftrC  •hown  1,1  figUrCS  2  and  3  1,1  whieh  ",e  45°  >>la,C 


(13 1 G)  corresponds  to  1  and  (he  G0°  plate  corresponds  to 


and  gives  the  following  relationship  between  A  and  B: 

'Pile  particular  solutions  for  uniform  loading  arc  given  by 

Vt'—ApX  i2 

<p2=BpxS 

where 

1  |  3a +1  ^ i _ 2  m  -^2_i  j 

t  A  "1— M  _ fi' 

2 8|772x?-,)x‘_'8xi— !)('!?‘2  v,Dc 

i _  ,, 

"»  4  |I„ 

8  (2X* —  1)X* — (SX2 —  1)(4 X* —  1 )  uc 


y 

-  | 

X«J 

x-'l 

3 

2.7031 

1  A6?l 

J.vn: 

a,«i7 

$  3S10 

4. 72oS 

-2.7031 

—1*  5*571 

-1.9137 

—10317 

-8.JU0 

-1.72v^ 

—.80(0  -.2*0  .saw  .1597  .01107  .MKB2  I 

-  01.557  -.02921  j  . .  .  -.07152  -.(MSI 

. 1 . o son  .00092 


Substitution  of  equations  (Bl7),  (1318),  and  (1319)  into  the 
equation 

tt>— yo+'/VT  +J/Vj 

gives  the  expression  for  tho  deflection  «>  of  the  plate  under  the 
assumption  of  parabolic  chordwise  deflection. 


DEFLECTION  AN'I)  STRESS  ANALYSIS  OK  THIN  SOLID  WINGS  OF  ARBITRARY  FLAN  FORM 
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DELTA  PLATE  OK  SYMMETRICAL  DOUBLE- WEDGE  AIRFOIL  SUCTION  WITH 
CONSTANT  THICKNESS  RATIO  UNDER  UNIFORM  LOAD 

For  r<  delta  plate  of  symmetrical  double-wedge  airfoil 
section  with  a  constant  thickness  ratio  the  thickness  is  a 
function  of  x  and  y  and  is  given  by  the  following  equations: 

/=  6<„f  (ogyg 

(xi-|)  (y!gi/gcx,) 

whore  /„  is  the  average  thickness.  From  these  expressions 
for  the  thickness  the  stiffness  can  bo  found  and  the  coolli- 
cicnts  in  tho  differential  equations  become 

27  Uc  , 

«i“— 


«i= 


27  Dc1 
*  8  '  " 


t )lJc*  . 

-  X,« 

o 


«<— 


Sl/Jc* 

SO 


-  x,' 


(1)20) 


2073  l)c> 

"s  ~  4480  -r‘ 

,  »75c*  , 

A|e3_so~  ■r' 

.  0/~c‘  . 
^  SO  x‘ 

,  201377c1  „ 

22400' 
yc"  ,, 

V  X' 


Solution  for  linear  chordwise  deflections. — lly  use  of  tho 
coeflieicnts  given  by  equations  (1120),  equations  (A22)  and 


(A  18)  for  linear  chordwise  deflections  may  bo  solved  for  <?t 
and  <i>q.  The  steps  in  the  solution  are  tho  same  in  form  ns 
those  for  tho  uniform-thickness  plato  and  tho  resulting 
equations  arc 


and 


40  /  l-x,T“! ,  .  \»/« 


(1121) 


1 

r  5 

20  7  2  j 

fxH-i  ) 

Y>o=27 

1  +'-x' 

y-f ' 

t  7~ 2  / 

/  °0  °\ 

/  -■ -q+^U-T'  log.x.+l-x,) 

V~4  7 


§  (1122) 


where 


~\'\2 


*)+80X* 


Solution  for  parabolic  chordwise  deflections. — By  use  of 
tho  coeflieicnts  given  by  equations  (B20),  equations  (A41), 
(A35),  and  (A33)  for  parabolic  chordwise  deflections  may  be 
solved  for  <pt,  yj*,  and  ip0.  The  stops  in  tho  solution  nro 
again  the  same  in  form  as  those  for  the  uniform-thickness 
plato  and  the  resulting  general  expressions  for  *>i,  yi,  and 
</>o  are 

c  T<_» 

Pi=2  log,Xi  (B23) 

7.-W 


Y2=  S  +  IK  —  (B24) 


C,xl'l'~2+C,xi  log,xrH7,x,-t-C,  (B25) 


where  tho  exponents  7,  arc  the  roots  of  the  characteristic  equation 

[I  (^"t) ~  -00x’ J=  0  (1120) 


For  n—  1,  2, . . .  G,  A„,  IJ„,  and  G„  aro  related  by 

«,=- 


or. 

y!-~-80y 

4 


ft— ■ 


f 

I  • 
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For  uniform  loud  tho  coefficients  in  llio  particular  integrals  of  equations  (B23),  (B24),  and  (B25)  are 


A  IQ  pi  ^+1024X^+320  x'-lG  (6+-^  X»)(10XM-1) 

/>c(20X*+l)(^-+1024Xl+320  X4)- 120(1  CX*+ 1)(1 0X‘+1) 

n  _40j£ _ _ ISdOX^D— 2(20X»  M)(54-^) _ 

2/  /Jc1(20x,+  ,)^8Vl  +  io24xl+32o|i^  X^-mOOX’+IKlOXH-l) 

n  _£  ,i _ §_  (.s^i  j.£  t*V  n  j.JL  PJl 

Lf  2 Ap  )5e  V  +2  1-)i/^+81  /} 

The  coefficients  At  to  Ae,  A„  C„  and  C,  arc  again  determined  bv  the  boundary  conditions  (A2G),  (ABC),  (A38),  and  (A39)  in 
which  tho  coefficients  given  by  equations  (B20)  nrc  substituted. 

For  Poisson’s  ratio  n  equal  to  l  and  X  =  -  =  1 ,  tbc  solution  of  the  characteristic  equation  (B20)  leads  to  two  real  values  and 
o  C 

two  pairs  of  complex  conjugate  values  for  y.  Tho  identity 

V±I>3V  cos  (6  log,  z,)±ix“  sin  (b  log,*,) 

was  therefore  used  to  transform  the  terms  involving  the  complex  conjugate  values  into  real  form.  If  I  and  u  — .*»  the 

c  3 

solution  is 

S55=|^r0.00-t070*1’',,,-0.004363*1,“scos  (2.825  log, *,)+ 0.000893*, ,OT‘ sin  (2.825  log,*,) +  0.000294  -1 
tic  L  ■r> J 

£-0.003890*,' ,,J  + 0.002 134r,,/1,s  cos  (2.825  log,  *,) -0.000381*,' mi  sin  (2.825  log,  *,)  +  0.01794  log,  *,  i  0.0017a'iJ 

v"o=Sj  j^0.000771 5*,‘-9,J- 0.0000708*, 1 cos  (2.825  log,*,)  +0;001234x,M-0B sin  (2.825  log, *,)  + 

0.03331*,  log,  *, -0.04090*, +0.04020 J 


APPENDIX  C 

NUMERICAL  PROCEDURE  FOR  SOLVING  DIFFERENTIAL  EQUATIONS 


In  oases  where  the  equations  of  the  present  theory  cannot 
be  solved  exactly  a  numerical  method  must  be  used.  In 
this  appendix,  equations  (Alt))  and  equations  (A34)  and 
(A35)  arc  set  up  in  difference  form  for  numerical  solution. 
Initially  the  assumption  is  matle  that  the  functions  involved 
in  the  differential  equations  are  continuous  and  nonsingular. 
In  this  case,  first  and  second  derivatives  can  be  expressed  by 
the  standard  difference  forms 

/< l‘y\ 

W7,~  V 

di a  _y.-4-y.-1 

dx)„  t 

whoro  t  is  the  distance  between  equally  spaced  station  points. 

In  the  following  development  five  equally  spaced  span- 
wise  stations  arc  used;  however,  the  extension  to  a  different 
number  of  stations  can  bo  readily  made. 


First,  consider  equation  (A19;  resulting  from  the  linear 
theory 

(hmT-  2(1  -p)(«,^,0,=7>1-(^££pi(x)  <fr,)"=7, 

Because  of  the  nature  of  the  tip  boundary  conditions  for 
this  equation,  it  can  be  conveniently  put  in  the  form 

?"= r„  (01) 

where 

7'=(6,v.T-2  (l-a)«iv.' 

In  finding  the  difference  cquution  equivalent  to  equation 
(Cl),  the  quantity  is  found  in  matrix  form,  from 

this  expression  is  subtracted  the  matrix  equivalent  of 
2(1—  ii)aiv>i',  the  resulting  expression  for  if’ is  multiplied  by 
a  differentiating matrix;  and  the  product  is  equated  to  the 
right-hand  side. 


The  quantity  (Aiv’i")'  at  the  half-stations  can  be  expressed  in  matrix  form  as  follows: 


Vl"o 

'1  -2  1 

Vi*l 

1  -2  1 

Vi  "2 

1 

1  -2  1 

“’i  ”3 

1  -2  1 

Vi". 

1-2  1. 

Vi-i 

Vio 

Vi  i 

V|2 

V|3 

Vn 

Vl5 


(C2) 


where  the  second  subscript  denotes  the  station  point,  the  subscript  at  the  root  station  being  0  and  at  the  tip  5.  The 
root  boundary  conditions  nrc  now  applied;  namely, 


Vi(0) — 0— -Vio 


V  i  1.0) =0—-- — t,( — 

Thus,  after  the  values  of  vj>o=0  and  vi-,~  Vn  arc  substituted,  equation  (C2)  becomes 


V."o 

-  2 

Vn 

Vi"i 

-2  1 

Vl2 

Vi  "1 

1 

1  -2  1 

Vn 

Vi  "i 

1  ~2  1 

Vn 

Vi"( 

1  -2  1: 

Vn 

•  ■ 


t. 

!  '• 


t  .  ■ 

t 

i 

\  • 
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(6iv»i")o 

(W). 

(61 

(4.V."). 

Ono  of  tho  tip  boundary  conditions  is 
'Inis, 

(4i»0!ia 

(4iv>i")«2 


r  o 
1—2  1 
1  -2  1 
1  -2  1 
1  -2 


(4lV>l*)i- 1=0— (4|(Pl")j 


(b^'%,,  =7 

(6.P.");/. 

(4nPi*)wj 

The  matrix  equivalent  for  tile  second  term  of  2’  is 
(uiVi'hrt  ■ 

(om')w  a 

20-m)  («,v»/)5«  =2Clrz£> 

(oifi')>n  L 

Therefore  T  becomes 


'-1  1 
—1  1 
-1  1 


'  (&.¥>.")  0 

1  (4iiPiff)a 

—  1  1  (tis?i")a 

-l-J  (&.*.'). 


1 

-1  1 
-1  1 
-1  1 


-1  1 _ I  Iv’isl 


Tn 

T,n 

Tin  =  ~ 

Tln 

Tin 


-1  1 
-1  1 
-1  1 
-1  1 
-1J 


2(1 -a)  I 


ir  2 
-2  1 
1  -2  1 
b»  1  -2 

4.JL  1 

’  1 
-1  1 
-1  1 
-1  1 
-1  1. 


The  right-hand  side  of  equation  (Cl)  can  now  be  equated  to  the  derivative  of  equation  (C5);  thus, 


'Ill 

'-11  1  Tm 

ll2 

_1 

—  1  1  Tm 

2u 

€ 

-1  1  Tin 

2h 

-1  lj  Tm 

Tin 

In  order  to  obtain  qa,  tho  boundary  condition 


at  x—l  must  be  used.  In  other  words,  T  goes  from  T,n  at  station  4^  to  0  at  station  5.  A  straight  line  drawn  between 
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o  Ttn 

these  two  points  would  havo  the  slope  —  The  value  of  qu  is  considered  to  be  this  slope;  therefore, 


■-1  1 

Tin 

-1  1 

Tin 

-1  1 

Tm 

t 

~1  1 

Tin 

-1. 

Tm 

{.»  =7. 
1u 

{is/2 


*1-2  1 
1  -2  1 
1  -2  1 
1  -2 
lj 


2 

-2  1 
1  -2  1 
1  -2  I 

1  -2  1. 


2(1  -a) 

; 


•-1  1 

-1  1 

-1  1 

-1  1 

-1. 


-1  1 

-1  1 

-1  1 


If  the  matrix  multiplication  is  carried  out,  tho  difference  equivalent  of  equation  (Cl)  finally  becomes 


hi 

¥>n 

•In 

={->  (0,1-2^  [£>,lj 

*>I2 

hi 

*>li 

hi 

<Pll 

{.S/2 

fll 

UJU 


2&,o+46,i+4u  26u  612 

—26,i— 2tij  6ii+46|2+6ii  — 26ij— 26u  6,2 

|Cj)=  in  — 2612 — 26u  4i2+46d+6h  — 26|j— 2bu 

in  — 26|j— 2611  6i2+4tu  —26k 

in  —  26m  &n  . 

•_aM/2--<*l,J/2  “l,J/2 

®I.3«  —  al,2/2— “l,i/2  °I.5A 

[Oil=  Ui,j«  ~~h,sn — di.m  “1,7/2 

_“l,7/2  —  r7!.  7/2  —  W  1,9/2  “l.9/2 

_  “l.9/3  — ai.9/r_ 

In  order  to  determine  y0  from  yj,,  use  must  bo  made  of  equation  (A18) 

or,  by  use  of  the  boundary  condition  n(0)*=p»'(0)=0, 


> 
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In  matrix  form  equation  (C7)  becomes 


Poi 

*1 

-1/2 

"l/"io 

‘1111  1' 

-11111” 

l»u 

Poi 

1  1 

1/2  l 

l/«n 

1111 

till 

/'I2 

pm 

1  1  1 

1/2  1  1 

1/«12 

1  1  1 

1  1  1 

/'|3 

Pm 

till 

1/2  1  1  1 

!/"■» 

1  1 

I  1 

Pu 

Poj 

.11111. 

.1/2  1111. 

1. 

/*n/2 

-1 

-1/2 

‘("j/",)o 

”  2 

P11 

1  1 

1/2  1 

(02/01)1 

-2  1 

P12 

1  1  1 

1/2  1  1 

(«2/«l)2 

1  -2  1 

Pis 

1111 

1/2  I  1  1 

("2/"l)3 

1  -2  1 

Pit 

.11111. 

.1/2  1111. 

(«2/«l)l. 

1  -2  1. 

Pis 

(CS) 


Tims,  if  the  values  of  q,  (which  can  he  determined  numerical!}  or  analytically  according  to  preference  and  feasibility;  are 
known,  the  values  of  pi  can  lie  found  by  solving  equation  (CO)  and  the  values  of  p0  in  turn  by  means  of  equation  (C8). 

The  foregoing  development  applies  to  the  case  where  only  linear  ehordwise  deformations  are  allowed.  A  similar  pro¬ 
cedure  is  followed  in  expressing  the  differential  equations  pertaining  to  the  parabolic  theory  in  difference  form,  only  the 
results  are  shown  herein. 


The  matrix  equivalent  to  equations  (A34)  and  (A35)  is 


where 


'/ii 

<Z.2 

'/l3 

I(<?1l_2(kriO|D), 

!» 

7is/2 

'/21 

f/22 

223 

~  [c,i-H‘T7^(z?2i 

plfti  8(17)[/;,i+ 

2« 

4(1— «*)I/iJ 

<Z«/2 

_ 

V=II 
P12 
Pi; 
Pi  i 
Pis 
P/i 
P  n 
Pa 
Pa 
P» 


r26.0+46,i+6,}  — 26„i—26,j  6,8 


-26„-26„ 

6,8 


6,1+46,8+6,3  —26,8—26,j 

“26,8—26,3  6,8+46,3+6,1 

6«  —26,3—26,, 
6,i 


^.3 

26,3  26,i  6„ 

6,i+46„  —26,( 

26,i  6,i  . 


(«,!= 


— ",,i/2“ "»,s/i 

",.3/3 


",.3  n 

— ",,3/2“",.3rj  ",.5/2 

",,5/2  “  ",,3«  —  ",.7/2  ",,7/2 

",,7/2  “  ",.7/8  —  ",.9/2 

",.9  // 


",,9/2 


“",.9/2.1 


(C9) 
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!£]= 


«H 


flu 


«|J 


and  Qt  and  q,  are  the  right-hand' sides  of  equations  (A34)  and  (A35),  respectively;  that  is, 

— erjw 

With  vi  and  v>i  known,  vo  can  be  obtained  by  use  of  equation  (A33) 

l/u io 


Voi 

Wj 

Voo 

Vm 

IvXm 


1 

1  1 

111,. 

1111 


1/2 

1/2 

1/2 

1/2 


1 

1  1 
1  l  1 


Li  i  i  i  iJLi/2  i  i  i 

(«»/«l)o 

(a, la,), 

(a, /a,), 

(«j/«i)s 


l/«n 


(<ij/«i).J 


1/Uii 


-2  1 
1  -2  1 
1  -2 
1 


1 /(in 


l/«iJ 


1111 
1  1  1 


r 

i 

i  i 
i  i 
i 


iiii 
i  i  i 
i  i 
i  i 

i 


"1 

2'n 

Pu 

Pi) 

Pu 

J 

Pul- 

2  1J 


W«.)o 


(aj/«i)i 


(«i/«i)i 


(a, la,). 


(a,la,)J 


■  o 

-2  1 
1-2  1 
1  -2  1 
1  -2 


(CIO) 


It  should  be  noted  that',  as  can  be  expected,  the  matrix  I 
equations  (CO)  and  (C8),  are  merely  special  cases  of  equations 
(G'9)  and  (CIO),  respectively.  In  addition,  the  square 
matrices  m  equations  (CO)  and  (C9)  are  symmetric,  a  result 
that  is  consistent  with  the  fact  that  the  (inferential  equations 
under  consideration  are  self-adjoint. 

In  the  beginning  of  tins  appendix  the  assumption  was  made 
that. the  functions  involved  in  the  differential  equations  are 
continuous  and  nonsingular.  The  difference  solution,  how¬ 
ever,  may  be  adequate  for  some  cases  in  which  this  assump¬ 
tion  is  not  strictly  correct.  For  instance,  the  deflections  of 
a  plate  with  a  discontinuous  stiffness  distribution  could  con¬ 
ceivably  bo  not  very  different  from  tho  deflections  of  a  plate 


I  w itli  a  continuous  stiffness  distribution  closely  approximating 
the  discontinuous  distribution  except  in  the  neighborhood  of 
tho  discontinuity.  The  results  yielded  by  the  difference  solu¬ 
tion  in  thisenso  would  be  those  associated  with  the  continuous 
stiffness  distribution.  The  number  of  stations  may  have  to 
be  increased,  however,  in  order  to  minimize  the  inuecuraey 
introduced  by  the  discontinuity  or,  in  other  eases,  by  a 
singularity.  The  case  of  the  symmetrical  double-wedge  air¬ 
foil  section,  coiistant-thickness-ratio  delta  plate,  discussed 
in  the  body  of  this  report,  is  an  example  of  a  treatment  of  a 
singularity.  In  this  case,  although  the  solution  is  singular, 
adequate  accuracy  is  obtained  by  the  difference  solution  if 
ten  equal  intervals  arc  used. 


■»- 


APPENDIX  D 

DEFLECTION  AND  STRESS  EXPERIMENTS  ON  SOME  TRIANGULAR  CANTILEVER  PLATES 


Test  specimen*. — Tlio  specimens  tested  were:  (1)  a  -15° 
right-triangular  plate  clamped  along  ono  leg  and  (2)  a  00° 
right-triangular  plate  clamped  along  the  longer  leg.  Each 
specimen,  cut  from  24S-T4  aluminum-alloy  sheet  of  0.250- 
inch  thickness,  had  a  length  perpendicular  to  the  clamped 
edgo  of  30  inches. 

Method  of  testing. — Figure  12,  a  photograph  of  the  test 
setup,  shows  the  methods  of  clamping,  loading,  and  measure¬ 
ment  of  deflections.  A  1,000,000-pound  clamping  load  (held 
constant  during  the  test)  was  applied  to  the  root  area  of  each 
specimen  and  a  uniform  load  of  0.204  psi  wan  applied  by 
2-ineh  washers  giving  a  tip  deflection  in  each  cane  of  approx¬ 
imately  •}  inch. 

The  deflections  were  measured  by  dial  gages  placed  at  the 
points  indicated  in  figures  2  and  3. 

Stresses  were  obtained  from  the  45°  specimen  only.  On 
this  specimen,  13  resistance  Wire  rosette  strait,  gages  were 
placed  at  the  points  indicated  in  figure  5.  Tlio  plate  was 
loaded  with  2-inch  washers  in  four  increments  of  0.0847  psi 
per  increment  and  tho  maximum  tip  deflection  was  1.13 
inches.  Headings  of  all  tho  strain  gages  were  recorded  at 
each  increment  of  loading. 

Analysis  and  discussion  of  data. — .The  deflection  to  was 
plotted  figures  2  and  3  in  terms  of  tho  nondimcnsional 


parameter  to  Dipl',  in  which  the  elastic  constants  were  taken 
as  A’=  10.GX  10*  psi  and  It  was  found  that  tho  dial- 

u 

gage  forces  reduced  the  tip  deflection  of  tho  plate  by  approx¬ 
imately  2  percent;  however,  since  this  error  is  of  tho  same 
order  of  magnitude  ns  thnt  in  the  mnlerial  properties  and 
from  other  sources,  no  corrections  are  mado  in  tho  results 
presented. 

The  readings  of  each  of  the  30  individual  strain  gages  were 
plotted  against  load,  and  the  slope  of  each  of  tho  resulting 
linear  curves  wns  taken  as  the  average  strain  per  unit  load 
of  tho  individual  gage.  The  principal  stresses  were  then 
calculated  nml  plotted  in  figure  5  in  terms  of  tho  nomlimen- 
sional  parameter  at ,lpl*. 
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